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Serum  penicillin  G  falls  to  low  levels  2  weeks  after  injection  as  benzathine  penicillin  G  (BPG)  in  young  adults.  Using  Pmetrics 
and  previously  reported  penicillin  G  pharmacokinetic  data  after  1.2  million  units  were  given  as  BPG  to  329  male  military  re¬ 
cruits,  here  we  develop  the  first  reported  population  pharmacokinetic  model  of  penicillin  G  after  BPG  injection.  We  simulated 
time-concentration  profiles  over  a  broad  range  of  pediatric  and  adult  weights  after  alternative  doses  and  dose  frequencies  to  pre¬ 
dict  the  probability  of  maintaining  serum  penicillin  G  concentrations  of  >0.02  mg/liter,  a  proposed  protective  threshold  against 
group  A  Streptococcus  pyogenes  (GAS).  The  final  population  model  included  linear  absorption  into  a  central  compartment,  dis¬ 
tribution  to  and  from  a  peripheral  compartment,  and  linear  elimination  from  the  central  compartment,  with  allometrically 
scaled  volumes  and  rate  constants.  With  1 .2  million  units  of  BPG  given  intramuscularly  every  4  weeks  in  four  total  doses,  only 
23.2%  of  5,000  simulated  patients  maintained  serum  penicillin  G  trough  concentrations  of  >0.02  mg/liter  4  weeks  after  the  last 
dose.  When  the  doses  were  1.8  million  units  and  2.4  million  units,  the  percentages  were  30.2%  and  40.7%,  respectively.  With  re¬ 
peated  dosing  of  1.2  million  units  every  3  weeks  and  every  2  weeks  for  4  doses,  the  percentages  of  simulated  patients  with  a  peni¬ 
cillin  G  trough  concentration  of  >0.02  mg/liter  were  37.8%  and  65.2%,  respectively.  Our  simulations  support  recommendations 
for  more  frequent  rather  than  higher  BPG  doses  to  prevent  recurrent  rheumatic  heart  disease  in  areas  of  high  GAS  prevalence  or 
during  outbreaks. 


Studies  performed  during  the  60  years  since  the  original  inves¬ 
tigations  (1)  of  the  persistence  of  serum  penicillin  G  after  par¬ 
enteral  administration  of  benzathine  penicillin  G  (BPG)  generally 
found  that  serum  concentrations  remain  above  a  putative  mini¬ 
mum  protective  level  for  the  prevention  and  treatment  of  group  A 
Streptococcus  pyogenes  (GAS)  infections.  This  minimum  threshold 
differs  among  authorities  but  is  usually  set  as  >0.01  to  0.03  p,g/ml 
(2)  for  between  3  and  4  weeks  after  a  dose. 

Recently,  we  reported  that  the  serum  penicillin  G  levels  in  329 
military  trainees  fell  more  rapidly  than  expected  during  the  29 
days  following  a  1.2-million-unit  intramuscular  injection  of  the 
only  formulation  of  BPG  currently  available  in  the  United  States 
(2).  That  study  was  prompted  by  observations  of  unexpectedly 
high  GAS  treatment  failures  with  BPG,  defined  as  a  failure  to  erad¬ 
icate  GAS  from  the  throat  (3-6),  and  previous  reports  of  unde¬ 
tectable  serum  penicillin  G  3  weeks  after  the  same  dose  (5,  7). 

In  this  report,  we  develop  a  population  pharmacokinetic  (PK) 
model  of  penicillin  G  concentrations  following  intramuscular  ad¬ 
ministration  of  1.2  million  units  of  BPG  in  healthy  young  adult 
males,  using  data  from  our  prior  study.  We  set  three  objectives:  (i) 
to  find  the  optimal  body  size  metric,  i.e.,  weight,  body  surface  area 
(BSA),  or  allometry,  for  scaling  PK  parameter  values;  (ii)  to  use 
the  final  model  for  Monte  Carlo  simulations  to  explore  alternative 
BPG  dosing  regimens  and  their  probabilities  of  maintaining  se¬ 
rum  penicillin  G  concentrations  of  >0.02  ptg/ml  (as  the  middle  of 
the  typical  target  range)  at  all  times  during  the  month  following  an 
injection  of  1.2  million  units;  and  (iii)  to  understand  the  kinetic 
distribution  of  penicillin  outside  the  serum  (i.e.,  in  tissues)  in 


relation  to  the  serum  time-concentration  profile.  These  objectives 
may  be  relevant  to  future  dosing  for  adults  and  children  for  ther¬ 
apy  of  both  GAS  upper  respiratory  tract  infections  and  secondary 
rheumatic  fever  prophylaxis  as  well  as  for  treatment  of  syphilis. 

MATERIALS  AND  METHODS 

Subjects.  As  described  previously  (2),  we  enrolled  329  male  military  train¬ 
ees  in  two  cohorts  and  monitored  them  for  29  days  after  their  enrollment 
and  injection  of  1.2  million  units  of  BPG.  None  of  the  subjects  were  aller¬ 
gic  to  penicillin. 

Procedures.  We  enrolled  the  first  cohort  in  January  2007  (165  subjects) 
and  the  second  in  February  2007  (164  subjects).  The  Naval  Health  Research 
Center  institutional  review  board  approved  the  study  (protocol  no. 
NHRC.2007.0022),  and  all  subjects  provided  written  informed  consent.  On 
day  0  for  each  subject,  we  gave  a  gluteal  intramuscular  dose  of  900  mg  of 
penicillin  G  as  2  ml  of  a  solution  containing  1 .2  million  units  of  BPG  (L-A 
BiriUin,  catalog  no.  1138883,  National  Drug  Code  [NDC]  60793070110; 
Monarch  Pharmaceuticals).  This  was  the  only  BPG  dose  that  each  subject 
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TABLE  1  Sampling  schedule  for  samples  2  and  3“ 


Group  ( n  =  12-18) 

Day  for  sample  2 

Day  for  sample  3 

i 

3 

16 

2 

4 

17 

3 

6 

20 

4 

7 

21 

5 

8 

22 

6 

9 

23 

7 

10 

24 

8 

13 

27 

9 

14 

28 

10 

15 

29 

a  All  subjects  were  dosed  on  day  0  and  sampled  on  day  1. 


received  in  the  study.  On  day  1,  which  was  sample  1,  we  withdrew  one  blood 
sample  from  all  subjects  by  venipuncture.  Samples  2  and  3  were  taken  on 
different  days  for  cohorts  of  12  to  18  subjects,  to  cover  the  entire  time-con¬ 
centration  profile,  as  shown  in  Table  1. 

Laboratory  analysis.  Samples  were  stored  at  — 70°C  for  up  to  18 
months  until  analysis.  These  samples  were  analyzed  in  batches  of  20  in 
random  order  to  remove  any  confounding  due  to  penicillin  G  degrada¬ 
tion.  We  used  our  previously  described  and  validated  liquid  chromatog¬ 
raphy/tandem  mass  spectrometry  (LC/MS/MS)  assay  (2). 

Population  PK  model  building.  For  population  PK  modeling,  we 
used  the  Non-Parametric  Adaptive  Grid  algorithm  (8)  in  the  Pmetrics 
package  for  R  (version  1.2.7)  (9)  to  fit  candidate  PK  models  to  the  time- 
concentration  data  for  penicillin  G.  For  all  models,  we  set  the  time  units  to 
be  days  rather  than  the  more  typically  used  hours. 

We  began  with  the  simplest  structural  model,  which  was  a  two-com¬ 
partment  model  that  had  a  depot/dosing  compartment  and  a  central  com¬ 
partment.  We  then  tested  a  model  with  an  additional  peripheral  tissue 
compartment.  Although  only  three  samples  were  obtained  from  each  sub¬ 
ject,  the  distribution  of  sampling  times  over  the  entire  time-concentration 
profile  allowed  us  to  at  least  test  these  two  models.  All  models  had  linear 
rate  constants  of  absorption  (KJ  from  dosing  to  the  central  compartment 
with  a  volume  of  distribution  ( V) ,  a  linear  elimination  rate  constant  from 
the  central  compartment  (Ke),  and,  when  applicable,  linear  transfer  rate 
constants  to  and  from  the  peripheral  compartment  (Kcp  and  KPC,  respec¬ 
tively).  We  assumed  that  measured  penicillin  G  concentrations  repre¬ 
sented  the  central  compartment.  Finally,  we  tested  the  scaling  of  PK  pa¬ 
rameters  to  three  different  size  models:  weight  ( V  =  V0  ■  wt),  body  surface 
area  [V  =  V0  ■  BSA),  or  allometric  size  (V  =  V0  ■  wt,  Ks  =  Ke0  ■  wt~0'25) 
(where  wt  is  weight  in  kg  normalized  to  70  kg  and  BSA  is  the  Mosteller 
body  surface  area  [10]  calculated  from  the  height  and  weight  of  each 
subject,  normalized  to  the  population  mean  of  1.94  m2). 

In  Pmetrics,  a  portion  of  the  random  error  in  the  model  is  attributed  to 
the  drug  assay,  with  a  multiplicative  error  term,  y ,  to  capture  process  noise 
associated  with  measurement  of  the  serum  concentration,  such  as  uncer¬ 
tainty  in  dose  times.  Pmetrics  weights  each  observation  by  the  reciprocal 
of  this  total  assay  variance  (Fisher  information),  calculated  as  (y  ■  SD)2, 
where  SD  =  C0  +  Q  •  [PCN]  +  C2  ■  [PCN]2  +  C3  ■  [PCN]3  (where  SD  is 
the  standard  deviation,  C0  to  C3  are  coefficients,  and  [PCN ]  is  the  mea¬ 
sured  penicillin  G  concentration).  We  fixed  C0  at  0.00025  as  half  the  assay 
limit  of  quantification  of 0.005  ng/ml,  Cj  at  0. 10,  and  C2  and  C3  at  0  for  an 
assay  with  a  10%  coefficient  of  variation. 

We  selected  the  final  model  on  the  basis  of  minimizing  the  Akaike 
information  criterion  (AIC)  (11),  which  is  a  function  of  the  likelihood  of 
the  model,  penalized  by  the  number  of  parameters  in  the  model.  Addi¬ 
tionally,  we  factored  bias  (mean  weighted  predicted-observed  error)  and 
imprecision  (bias-adjusted,  mean  weighted  squared  predicted-observed 
error)  into  the  selection  of  the  final  model. 

Visual  predictive  checks.  To  provide  some  assessment  of  the  ability  of 
the  final  model  to  accurately  represent  the  study  population,  we  used  the 
simulation  module  of  Pmetrics  to  create  5,000  sets  of  PK  parameter  values 


TABLE  2  Characteristics  of  study  subjects  (n  =  329) 


Variable 

Mean  value  (SD) 

Range 

Age  (yr) 

20  (1.9) 

17-32 

ht  (cm) 

176.8  (7.8) 

144.8-195.6 

wt (kg) 

76.7(11.6) 

50.0-109.1 

BSA  (m2) 

1.94  (0.17) 

1.53-2.42 

by  randomly  selecting  them  from  the  probability-weighted  normalized 
distribution  and  covariance  of  values  for  each  parameter  in  the  model,  a 
technique  known  as  Monte  Carlo  simulation  (12).  As  in  the  study,  each 
simulated  adult  was  administered  1.2  million  units  of  BPG  once  intra¬ 
muscularly,  and  penicillin  G  concentrations  were  simulated  daily  there¬ 
after  for  1  month.  We  calculated  and  plotted  the  10th,  50th,  and  90th 
percentiles  of  the  concentrations  of  penicillin  G  in  all  5,000  simulated 
patients  versus  time.  Superimposed  upon  these  plots  were  the  actual  con¬ 
centrations  of  penicillin  G  measured  in  the  329  actual  subjects.  This  visual 
predictive  check  (13)  was  considered  good  if  the  distribution  of  concen¬ 
trations  in  the  simulated  population  was  similar  to  that  in  the  actual 
population. 

Probability  of  target  attainment.  Using  the  final  model,  the  simula¬ 
tor,  and  the  probability  of  target  attainment  (PTA)  functions  in  Pmetrics, 
we  explored  various  alternative  dosing  regimens  for  BPG.  For  each  dosing 
regimen,  we  simulated  5,000  patients.  Pmetrics  allows  simulation  with 
covariates,  using  the  correlation  matrix  of  Bayesian  posterior  parameter 
values  and  covariates  in  the  study  population,  together  with  the  specified 
mean  and  standard  deviation  of  any  co  variates,  all  to  calculate  the  correct 
covariate-parameter  covariance  matrix  for  simulations.  In  this  case,  the 
only  covariate  was  weight,  which  we  set  for  simulations  to  a  mean  of  50  kg 
and  a  standard  deviation  of  50  kg,  truncating  simulated  values  to  be  within 
the  range  of  25  to  1 10  kg.  This  resulted  in  a  uniform  distribution  of  weight 
between  our  specified  limits,  i.e.,  a  homogeneously  mixed  group  of  sizes 
corresponding  to  children,  adolescents,  and  adults.  We  used  the  popula¬ 
tion  distributions  for  the  model  pharmacokinetic  parameters,  limiting 
simulated  values  to  between  0  and  5  times  the  maximum  in  order  to 
simulate  a  more  diverse  population  than  the  original  study  population. 
For  each  set  of  simulated  pharmacokinetic  parameters  and  weights,  we 
calculated  the  proportion  of  simulated  patients  who  maintained  a  serum 
penicillin  G  concentration  of  >0.02  mg/liter  at  all  times  during  the  sub¬ 
sequent  28  days.  For  doses  of  >1.2  million  units,  we  assumed  that  con¬ 
centrations  would  increase  proportionally  to  the  dose,  since  we  did  not 
have  data  for  higher  doses. 

Finally,  since  the  final  model  included  a  peripheral  tissue  compart¬ 
ment,  we  sought  to  understand  how  much  penicillin  G  kinetically  distrib¬ 
utes  from  the  serum  to  extravascular  tissue  and  if  penicillin  G  accumulates 
in  bodily  fluids  outside  the  serum  over  time  with  repeated  BPG  dosing.  To 
do  this,  we  simulated  one  additional  population  of  5,000  adults  who  re¬ 
ceived  1.2  million  units  of  BPG  monthly  for  6  months,  to  achieve  steady 
state.  From  this  population,  we  simulated  the  amount  of  penicillin  G  in 
the  peripheral  compartment  at  two  times:  1  month  after  the  first  dose  and 
again  after  the  final  (sixth)  dose.  For  each  of  the  5,000  simulated  adults  at 
each  of  the  two  time  points,  we  calculated  the  relative  distribution  of 
penicillin  G  to  the  peripheral  compartment  as  the  ratio  of  the  model- 
calculated  amount  of  penicillin  G  in  the  peripheral  compartment  to  the 
amount  in  the  central  (serum)  compartment.  We  also  assessed  accumu¬ 
lation  of  penicillin  G  in  the  peripheral  compartment  over  time  as  the  ratio 
of  the  amount  of  penicillin  G  1  month  after  the  final  dose  to  the  amount  1 
month  after  the  first  dose. 

RESULTS 

Characteristics  of  the  study  population  are  shown  in  Table  2.  The 
distribution  of  penicillin  G  concentrations  per  day  is  included  in 
Fig.  1  (circles).  There  was  one  outlier  subject  whose  measured 
penicillin  G  concentration  on  day  22  was  0.14  (xg/ml,  versus  a 
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FIG  1  Distribution  of  measured  penicillin  G  concentrations  (circles)  in  the 
population.  Lines  are  the  indicated  percentiles  of  5,000  simulated  concentra¬ 
tion-time  profiles.  The  gray  shading  around  the  percentile  lines  represents  the 
95%  confidence  interval  around  each  percentile.  The  dotted  horizontal  line  at 
0.02  mg/liter  is  the  suggested  minimum  protective  concentration  of  penicillin 
G  against  group  A  streptococcus.  Note  that  the  majority  of  measured  concen¬ 
trations  fall  below  this  threshold.  The  dashed  horizontal  line  is  the  limit  of  the 
assay,  which  is  below  all  measured  concentrations  but  is  above  some  simulated 
values  at  the  end  of  the  dosing  interval.  As  a  visual  predictive  check  (13)  of  the 
model,  the  distribution  of  the  simulated  profiles  is  similar  to  that  of  the  ob¬ 
served  concentrations,  suggesting  that  the  model  describes  the  data  well. 


mean  of  0.006  |ag/ml  for  the  other  subjects.  However,  we  retained 
this  outlier  measurement  in  the  analysis,  as  we  did  not  have  a 
defined  reason  to  exclude  it. 

Population  model.  We  compared  AIC,  bias,  and  imprecision 
for  five  models  after  allowing  the  population-fitting  algorithm  to 
iterate  1,000  times  toward  the  convergent,  maximally  likely  distri¬ 
bution  of  parameter  values.  Details  of  the  models  and  their  AIC 
values  are  shown  in  Table  3.  On  the  basis  of  the  AIC,  bias,  impre¬ 
cision,  and  stronger  a  priori  justification  for  scaling  to  children,  we 
chose  model  5,  which  had  a  peripheral  tissue  compartment  and 
allometrically  scaled  central  compartment  volume  and  elimina¬ 
tion  from  that  compartment.  We  then  allowed  model  5  to  cycle 


TABLE  4  Population  parameter  values  after  a  single  dose  of  1 .2  million 
units  of  benzathine  penicillin  G  in  adults0 


Parameter 

Unit  of  measure 

Value 

Mean  (SD) 

Median 

K'O 

l/(day  •  70  kg0'25) 

4.11  (2.22) 

3.96 

K„ 

1/day 

0.40  (0.25) 

0.32 

Kcpo 

l/(day  •  70  kg0  25) 

6.91  (6.30) 

4.71 

^PCO 

l/(day  •  70  kg0  25) 

0.27  (0.34) 

0.12 

Vo 

Liters/70  kg 

260.96  (96.45) 

239.43 

Calculated  from  full- fitted 

concn  profiles 

AUC^ 

mg  ■  h/liter 

19.33  (6.09) 

18.68 

Clearance 

Liters/h/kg 

0.68  (0.24) 

0.65 

Half-life  of  absorption 

h 

72.55  (202.21) 

50.79 

Half-life  of  elimination 

h 

6.03  (4.29) 

4.19 

r 

uma x 

mg/liter 

0.14  (0.09) 

0.13 

T 

max 

h 

22.90  (32.10) 

9.60 

a  AUC0_oo,  area  under  the  concentration-time  curve  from  0  h  to  infinity;  Ka,  absorption 
from  dosing  to  the  central  compartment;  Ke ,  elimination  from  the  central 
compartment;  KCP ,  transfer  from  the  central  to  the  peripheral  compartment;  KPC> 
transfer  from  the  peripheral  to  the  central  compartment;  V,  volume  of  the  central 
compartment;  Cmax,  maximum  concentration;  Tmax,  time  to  maximum  concentration. 


until  convergence,  which  took  9,429  cycles.  The  final-cycle-esti- 
mated  y  from  the  pooled  participant  concentration  data  was  1.08. 
A  value  of  1.0  indicates  that  there  is  no  additional  process  noise  in 
the  study,  such  as  errors  in  the  recorded  times  of  sampling,  and 
that  the  study  procedures  were  carried  out  very  precisely. 

The  population  parameter  value  distributions  for  the  final 
model  are  summarized  in  Table  4,  and  the  full  marginal  distribu¬ 
tions  are  shown  in  Fig.  2.  The  parameter  values  appeared  to  most 
closely  approximate  a  log-normal  distribution.  The  half-life  of 
penicillin  G  absorption  after  dosing  as  BPG  is  much  longer  than 
the  half-life  of  elimination,  as  would  be  expected.  Using  the  pa¬ 
rameter-value  distributions  in  Table  4,  the  visual  predictive  check 
of  5,000  concentration-time  profiles  versus  the  observed  penicil¬ 
lin  G  concentrations  in  the  study  population  (329  subjects)  is 
shown  in  Fig.  1.  The  simulated  distribution  of  concentrations 
matches  the  observed  distribution  well,  suggesting  that  the  model 
describes  the  study  data  adequately  and  can  be  used  for  meaning¬ 
ful  PTA  analysis. 

Probability  of  target  attainment.  In  addition  to  the  standard 
dose  of  1.2  million  units  of  BPG  every  4  weeks,  we  studied  simu- 


TABLE  3  Model  statistics  after  1,000  cycles0 

Population  prediction  Individual  prediction 


Model 

Parameters 

AIC  value 

Bias 

Imprecision 

Bias 

Imprecision 

1 

K. „  K e,  V 

-2,362 

0.054 

6.15 

-0.39 

0.65 

2 

Ka>  Ke,  V,  Kcp,  Kpc 

-2,435 

4.81 

125.28 

-0.01 

1.18 

3 

Ks,Ke,V=  V0  ■  wt,  KCP,  Kpc 

-2,444 

4.98 

139.54 

-0.06 

0.97 

4 

Ks,Ke,V=  V0  ■  BSA,  KCP,  KPC 

-2,441 

5.05 

127.92 

0.08 

2.20 

5 

K„  K  =  KJ wt0  25,  V  =  V„  •  wt,  Kc P  =  iCCp/wt°-25, 

XPC  =  Ifpc/wt025 

-2,440 

3.51 

108.65 

-0.01 

1.16 

a  All  size-scaled  models  (models  3  to  5)  are  similar,  but  model  5,  with  allometric  scaling  for  body  size,  is  preferable  based  on  minimization  of  AIC  and  favorably  low  bias  and 
imprecision.  AIC,  Akaike  information  criterion,  with  the  lowest  value  indicating  the  most  likely  model;  bias,  mean  weighted  error  of  predictions  minus  observations;  BSA,  body 
surface  area  in  m2  normalized  to  a  mean  population  BSA  of  1.94  m2;  imprecision,  bias-adjusted  mean  weighted  squared  error  of  predictions  minus  observations;  Ka,  absorption 
from  dosing  to  the  central  compartment;  Kcp,  transfer  from  the  central  to  the  peripheral  compartment;  Ke,  elimination  from  the  central  compartment;  Kpc,  transfer  from  the 
peripheral  to  the  central  compartment;  V,  volume  of  the  central  compartment;  wt,  weight  in  kg  normalized  to  70  kg. 
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FIG  2  Marginal  distributions  ofparameter  values  in  the  final  model.  Ke0,  elimination  ofBPG  from  the  central  compartment  scaled  to  body  weight  (h  1  ■  kg  °'25); 
fCa,  absorption  from  muscle  to  the  central  compartment  (h  ]);  Kcp  and  Kpc,  distribution  ofBPG  from  the  central  compartment  to  the  peripheral  compartment 
and  back  (h_1);  V0,  volume  of  the  central  compartment  (liters/kg). 


lated  doses  of  1.8  million  or  2.4  million  units  every  4  weeks  as  well 
as  1.2  million  units  given  every  2  or  3  weeks.  The  results  are  shown 
in  Table  5  and  Fig.  3.  All  of  the  regimens  resulted  in  substantial 
proportions  of  the  study  population  with  concentrations  pre¬ 
dicted  to  fall  below  0.02  mg/liter  of  penicillin  G  in  the  serum  at  the 
end  of  the  dosing  interval,  with  only  23.2%  having  concentrations 
above  this  threshold  4  weeks  after  the  fourth  dose  of  1.2  million 
units.  The  regimen  with  the  highest  success  rate  was  1.2  million 
units  ofBPG  every  2  weeks,  with  65.2%  of  patients  having  con¬ 
centrations  of  >0.02  mg/liter  2  weeks  after  the  previous  dose  at 
steady  state. 


TABLE  5  Proportion  of  5,000  simulated  subjects  weighing  between  25 
and  110  kg  with  trough  serum  penicillin  G  concentrations  of  >0.02  mg/ 
liter  for  various  dosing  regimens 


Dose 

%  of  subjects  with  trough 
concnfl  of  >0.02  mg/liter 

1.2  million  units  monthly 

23.2 

1.2  million  units  every  3  wk 

37.8 

1.2  million  units  every  2  wk 

65.2 

1.8  million  units  monthly 

30.2 

2.4  million  units  monthly 

40.7 

a  The  trough  concentration  is  defined  as  the  serum  penicillin  G  concentration  just  prior 
to  the  next  dose. 


Accumulation  and  peripheral  compartment  distribution  of 
penicillin  G.  In  our  model,  4  weeks  after  the  first  dose  of  1.2 
million  units  ofBPG,  the  serum  penicillin  G  concentration  in  the 
simulated  adult  population  (using  a  median  population  weight  of 
74  kg)  was  a  median  of  0.005  mg/liter  (interquartile  range  [IQR], 
0.001  to  0.01  mg/liter).  For  comparison  as  a  validation,  in  the  real 
subjects,  this  value  was  very  similar:  0.004  mg/liter  (IQR,  0.002  to 
0.009  mg/liter).  In  the  peripheral  compartment  of  the  simulated 
population,  the  penicillin  G  concentration  was  0.02  mg/kg  of  body 
weight  (IQR,  0.005  to  0.04  mg/kg).  Note  that  this  does  not  corre¬ 
spond  to  an  actual  concentration  in  any  specific  tissue,  but  it  is  a 
weight-normalized  amount  of  drug  that  has  kinetically  distrib¬ 
uted  outside  the  measurable  serum;  i.e.,  this  is  a  mathematical 
phenomenon  used  to  explain  the  shape  of  the  observed  serum 
concentrations  with  respect  to  time.  Nevertheless,  we  examined 
the  characteristics  of  the  peripheral  distribution  of  penicillin  to 
determine  if  the  drug  was  likely  to  be  accumulating  outside  the 
serum.  The  median  ratio  of  the  peripheral  to  the  serum  penicillin 
G  concentration  after  4  weeks  was  3.8  (IQR,  2.6  to  5.1).  The  ratio 
of  the  peripheral  concentration  to  the  serum  concentration  4 
weeks  after  the  sixth  dose  was  3.7  (IQR,  2.7  to  5.0),  the  same  as  4 
weeks  after  the  first  dose,  which  indicates  that  the  drug  amount  in 
peripheral  tissues  is  not  larger  after  six  doses  than  after  only  one 
dose.  This  lack  of  substantial  accumulation  over  time  is  further 
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FIG  3  Proportions  of  simulated  patients  with  penicillin  G  trough  concentrations  of  >0.02  mg/liter  for  various  dosage  regimens  of  BPG  on  a  given  day  after  the 
fourth  dose,  i.e.,  at  steady  state. 


supported  by  comparing  the  median  ratio  of  the  serum  penicillin 
G  concentration  4  weeks  after  the  sixth  dose  to  the  serum  penicil¬ 
lin  G  concentration  4  weeks  after  the  first  dose,  which  was  1.23 
(IQR,  0.97  to  1.73).  In  the  peripheral  compartment,  this  ratio  was 
1.08  (IQR,  1.01  to  1.40).  In  other  words,  in  both  serum  and  pe¬ 
ripheral  tissues,  the  penicillin  G  concentration  4  weeks  after  the 
sixth  dose  is  similar  to  the  concentration  4  weeks  after  the  first 
dose. 

Together,  these  data  show  that  there  is  very  little  accumulation 
of  penicillin  G  in  the  serum  or  peripheral  tissue  compartment  over 
time  at  the  dose  used  in  this  study.  This  is  because  the  drug  is 
nearly  fully  cleared  from  the  body  with  monthly  dosing,  and  the 
weight-normalized  amount  of  penicillin  G  kinetically  distributed 
to  the  peripheral  compartment  is  stable  with  respect  to  time  at 
^2.5  to  5  times  the  serum  concentration. 

DISCUSSION 

We  present  the  first  reported  population  model  of  penicillin  G 
administered  as  BPG.  This  model  was  developed  with  >300 
young,  active,  and  healthy  adults  who  each  contributed  three  se¬ 
rum  samples  covering  the  period  from  the  first  day  after  adminis¬ 
tration  and  weekly  thereafter  until  1  month  after  dosing.  To  ex¬ 
plore  the  dose-exposure  relationships  with  the  currently  available 
preparation  of  BPG,  we  scaled  the  penicillin  G  PK  parameters  to 
weight,  BSA,  and  allometric  size  and  found  them  to  be  roughly 
equivalent.  In  our  model,  however,  allometric  scaling  was  slightly 
better,  and  it  allowed  for  better  extrapolation  from  adult  to  pedi¬ 
atric  populations  (14-16). 

It  is  clear  from  our  data  that  with  the  exception  of  biweekly 
BPG  dosing,  the  majority  of  patients  will  not  sustain  serum  peni¬ 


cillin  G  concentrations  above  the  MIC  of  >0.02  mg/liter  for  GAS 
during  the  entire  dosing  interval.  Here  we  have  demonstrated  this 
in  young  healthy  adults,  and  by  modeling  and  simulation,  we  ex¬ 
trapolated  our  conclusion  to  patients  with  weights  ranging  be¬ 
tween  25  and  110  kg.  This  extrapolation  appears  valid,  since  our 
predictions  of  the  proportions  of  patients  with  concentrations 
falling  below  0.02  mg/liter  3  to  4  weeks  after  doses  of  1.2  million  to 
2.4  million  units  of  BPG  (Table  4)  are  in  close  agreement  with 
observations  made  by  others  (3,  17). 

Penicillin  G  is  a  “time-dependent”  beta-lactam  antibiotic, 
which  is  a  class  of  antibiotics  whose  antibacterial  kill  rates  are 
maximized  when  the  serum  concentration  is  above  the  MIC  of  the 
organism  for  30  to  70%  of  the  dosing  interval  (18),  depending  on 
the  drug  and  the  organism.  However,  this  PK/pharmacodynamic 
(PD)  linkage  is  really  relevant  only  to  drugs  that  are  administered 
daily  or  multiple  times  a  day.  There  is  no  such  time-dependent 
model  of  the  effects  of  repository  penicillin  G  administered  as 
BPG  on  GAS.  Prolonged  concentrations  below  the  MIC  of  the 
organism  for  1  to  2  weeks  are  likely  not  optimal,  as  suggested  by  a 
study  over  a  32-week  period  from  November  1956  to  June  1957, 
where  the  breakthrough  rate  for  GAS  infections  in  military  re¬ 
cruits  within  the  month  after  administration  of  1.2  million  units 
of  BPG  averaged  3.7  per  100  patient-years  (19),  and  almost  all 
GAS  infections  occurred  >2  weeks  after  injection,  when  concen¬ 
trations  were  lower.  In  1996,  Lue  et  al.  reported  7.5  and  12.7 
breakthroughs  per  100  patient-years  among  children  and  adults 
receiving  1.2  million  units  of  BPG  every  3  weeks  and  every  4  weeks, 
respectively  (20);  i.e.,  those  with  longer  interdose  intervals  had 
more  breakthroughs  ( P  <  0.01).  In  a  comparison  of  twice- 
monthly  versus  once-monthly  BPG,  Kassem  et  al.  reported  3.7 
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breakthroughs  per  100  patient-years,  and  50%  of  the  infections 
occurred  during  the  third  or  fourth  week  after  injection,  versus 
only  21%  in  the  first  2  weeks  (21).  Despite  the  use  of  BPG  formu¬ 
lations  from  different  manufacturers,  which  can  affect  serum  pen¬ 
icillin  G  concentrations  (22) ,  these  studies  clearly  suggest  that  pro¬ 
tection  against  GAS  wanes  after  the  first  2  weeks  following  a  dose, 
which  is  consistent  with  our  own  observations. 

Waning  serum  penicillin  G  concentrations  may  also  be  re¬ 
flected  in  a  relatively  high  rate  of  failure  to  eradicate  GAS  from  the 
pharynx.  One  study  found  a  37%  rate  of  failure  to  eradicate  pha¬ 
ryngeal  GAS  10  to  14  days  or  29  to  31  days  after  a  single  injection 
of  600,000  units  (if  patients  weighed  <60  lb)  or  1.2  million  units 
(if  patients  weighed  >60  lb)  of  BPG  in  271  children  who  presented 
with  acute  pharyngitis  and  who  were  culture  positive  for  GAS  (6). 

Despite  this  evidence  that  waning  penicillin  G  concentrations 
are  linked  with  increased  numbers  of  breakthrough  infections  and 
failure  to  eradicate  GAS  in  those  patients  with  established  pharyn¬ 
geal/throat  infections,  the  overall  failure  rates  in  the  setting  of  BPG 
injections  monthly  or  every  3  weeks  are  still  relatively  low  com¬ 
pared  with  the  proportion  of  the  population  that  has  concentra¬ 
tions  below  the  MIC  threshold  of  0.02  mg/liter  by  2  weeks  after  an 
injection.  We  have  several  hypotheses  regarding  why  the  ability  of 
BPG  to  protect  against  GAS  infections  appears  to  be  better  than 
one  would  expect  from  the  associated  serum  penicillin  G  concen¬ 
trations.  First,  the  MIC  for  any  particular  GAS  isolate  may  be 
<0.02  mg/liter,  since  the  MIC  for  90%  of  >4,000  strains  was 
reported  to  be  <0.06  mg/liter  in  2004  (23).  Only  one  strain  out  of 
282  had  an  MIC  of  >0.012  mg/liter  (0.024  mg/liter)  in  1992  (24). 
No  isolates  had  an  MIC  of  >0.01  mg/liter  in  1965  (25).  More  than 
twice  as  many  patients  are  expected  to  maintain  serum  penicillin 
G  concentrations  of  >0.0 1  mg/liter,  compared  to  0.02  mg/liter  for 
1  month  (data  not  shown). 

The  second  hypothesis  is  that  the  young  healthy  military  males 
in  this  study  may  have  cleared  penicillin  G  faster  than  other  pop¬ 
ulations  after  injection  as  BPG.  Other  patients  may  retain  higher 
concentrations  for  longer  periods.  In  one  of  the  only  reported 
clinical  studies  of  the  same  doses  that  we  used  for  our  simulations 
(1.2  million,  1.8  million,  and  2.4  million  units),  higher  propor¬ 
tions  (10  to  20%)  of  adolescents  and  young  adults  maintained 
therapeutic  serum  penicillin  G  concentrations  2,  3,  and  4  weeks 
after  injection  than  the  subjects  in  our  study,  although  that  study 
used  a  bioassay  to  measure  penicillin  G  concentrations  (26). 

The  third  hypothesis  is  that  tissue  penicillin  G  concentrations 
are  likely  to  be  different  in  blood  (17,  27).  Since  most  GAS  infec¬ 
tions,  such  as  tonsillitis,  occur  in  tissues  other  than  blood,  we 
attempted  to  relate  the  serum  kinetics  of  penicillin  G  to  the  kinet¬ 
ics  outside  the  measurable  serum  compartment  over  time.  Specif¬ 
ically,  we  tried  to  address  the  question  of  whether  there  is  kinetic 
accumulation  of  penicillin  G  with  repeated  dosing  such  that  tissue 
concentrations  might  persist  after  the  clearance  of  drug  from  the 
serum,  extending  the  functional  efficacy  of  the  drug  to  prevent 
GAS  infection  and/or  disease.  However,  our  data  suggest  that 
there  is  little  residual  drug  in  the  serum  or  peripheral  compart¬ 
ments  by  the  end  of  a  28-day  dosing  interval,  even  with  monthly 
dosing;  that  is,  the  drug  is  nearly  fully  cleared  from  the  body 
within  1  month.  In  our  model,  we  did  find  that  the  average  con¬ 
centration  of  penicillin  G  in  the  peripheral  compartment  1  month 
after  a  dose  of  1.2  million  units  of  BPG,  when  normalized  to  body 
weight,  is  nearly  4-fold  higher  than  that  in  the  serum;  again,  how¬ 
ever,  this  is  a  kinetic  observation,  and  one  cannot  draw  specific 


conclusions  about  the  concentrations  of  penicillin  G  in  any  spe¬ 
cific  tissue.  For  example,  Peloso  et  al.  reported  that  the  average 
measured  concentration  ratio  of  tonsillar  to  serum  penicillin  G  in 
children  was  only  about  33%  for  the  first  2  weeks,  and  by  2 1  days, 
70%  of  children  did  not  have  detectable  tonsillar  penicillin  G  de¬ 
spite  having  measurable  serum  penicillin  G  concentrations  (17). 

The  fourth  hypothesis  is  that  GAS  disease  rates  depend  upon 
not  only  the  concentration  of  penicillin  G  in  the  serum  but  also  a 
combination  of  location-  and  season-dependent  GAS  prevalence, 
individual  infection  rates,  and  risk  of  progression  from  infection 
to  disease  (28-31).  Even  for  a  patient  with  subtherapeutic  serum 
penicillin  G  concentrations,  it  is  not  100%  certain  that  GAS  infec¬ 
tion  and  disease  will  follow.  However,  in  areas  with  higher  GAS 
prevalence,  more  frequent  BPG  injections  are  recommended, 
consistent  with  our  observations  of  the  waning  of  serum  penicillin 
G  levels  at  the  end  of  the  monthly  dose  interval  (4,  20,  32).  Future 
work  to  combine  our  PK  model  with  GAS  PK/PD  modeling  (33, 
34)  and  epidemiology  may  be  able  to  more  quantitatively  predict 
optimal  BPG  dosing  frequency  based  on  local  GAS  prevalence. 

In  summary,  we  have  developed  a  population  PK  model  to 
calculate  serum  penicillin  G  exposure  after  intramuscular  BPG 
dosing  in  healthy  young  adults.  Although  this  model  also  appears 
to  be  relevant  to  children,  there  is  a  paucity  of  actual  supportive 
pediatric  PK  data  (17).  With  the  current  dose  of  1.2  million  units, 
serum  concentrations  of  penicillin  G  will  fall  below  a  threshold  of 
0.02  mg/liter  in  the  majority  of  adults  and  likely  in  children  after  2 
weeks.  Higher  monthly  doses  seem  unlikely  to  substantially  in¬ 
crease  the  proportion  of  patients  whose  concentrations  remain 
above  0.02  mg/liter  for  the  entire  dosing  interval,  but  dosing  every 
3  weeks  or  even  every  2  weeks  will  be  more  successful.  However,  it 
may  well  be  that  a  minimum  level  of  0.02  mg/liter  is  unnecessarily 
high;  answering  this  question  will  require  further  mechanistic  and 
epidemiological  PK/PD  modeling  of  penicillin  G  and  GAS  disease. 
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